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ABSTRACT 

We present 15-GHz follow-up radio observations of eleven Swift gamma-ray burst 
(GRB) sources, obtained with the Arcminute Microkelvin Imager Large Array (AMI- 
LA). The initial follow-up observation for each source was made in a fully automated 
fashion; as a result four observations were initiated within five minutes of the GRB 
alert timestamp. These observations provide the first millijansky-level constraints on 
prolonged radio emission from GRBs within the first hour post-burst. While no radio 
emission within the first six hours after the GRB is detected in this preliminary analy- 
sis, radio afterglow is detected from one of the GRBs (GRB120326A) on a timescale of 
days. The observations were made as part of an ongoing programme to use AMI-LA as 
a systematic follow-up tool for transients at radio frequencies. In addition to the pre- 
liminary results, we explain how we have created an easily extensible automated follow- 
up system, describing new software tools developed for astronomical transient alert 
distribution, automatic requesting of target-of-opportunity observations, and robotic 
control of the observatory. 

Key words: gamma-ray burst: general - methods: miscellaneous - methods: obser- 
vational - virtual observatory tools - instrumentation: interferometers 



1 INTRODUCTION 



Since the discovery of GRB afterglow emission at X-ray 



(Costa et al. 19971, optical (van Paradijs et al. 19971 and 



radio frequencies (Frail et al. 19971, it has been demon- 



strated that broadband observations across the electromag- 
netic spectrum can be utilized to determine the micro- and 
macro-physical parameters of the GRB explosion and its im- 
mediate environment (e.g. Sari et al.|1998 Wijers & Galama 
1999). After the initial detection in gamma- and X-rays, 



rapid follow-up observations are often essential to capture 
the GRB afterglow in the optical, before it fades. To this 
end, automated, robotic optical observations from smaller 



telescopes are now routine (e.g. Akerlof et al. 2003 Woz- 



niak et al. 20051. In contrast, early follow-up observations 



of GRBs at radio wavelengths are rare (although some at- 
tempts have been made, see e.g. |Green et al. 1995; Koranyi| 



et al.|1995||Dessenne et al.|1996| 

In this paper we explain how we have achieved fully 
automated rapid follow-up observations of GRB events with 
the Arcminute Microkelvin Imager Large Array (AMI- LA), 



a 15-GHz aperture synthesis radio telescope located near 
Cambridge, UK. A primary beam width of 5.5 arcminutes 
and a sensitivity of ~ 3mJy s _1//2 make this an ideal facility 
for follow-up of the well localised GRB events observed by 



the Swift burst alert telescope (Swift-BAT, Gehrels et al. 
2004[|BartIielmy et al.|2005[ ). 



This paper is structured as follows: First we give the 
scientific motivation for obtaining early-time radio follow- 
up observations (Section A brief review of the relevant 
models for GRB emission is given, focusing upon predictions 
for radio emission at the time of the gamma ray emission and 
within the first few hours thereafter. We then explain how 
this programme will provide novel data, complementary to 
existing GRB follow-up catalogues. In Sections [3] - [5] we 
describe the AMI Large Array; the system we have put in 
place to initiate automated follow-up; and our data reduc- 
tion techniques. In Section [6] we report on system response 
times and present 15-GHz lightcurves. Finally, we describe 
our plans for further analysis of the AMI-LA data, and pos- 
sibilities for follow-up with alternative facilities, in Section 

m 
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2 SCIENTIFIC MOTIVATION 

2.1 GRB radio emission at short and long 
timescales 

The 'relativistic fireball' is the prevalent model to describe 
GRB emission across the electromagnetic spectrum. In this 
model a blast wave propelled by an ultra-relativistic jet 
ploughs through an external medium. Electrons accelerated 
by the shock emit synchrotron radiation in the magnetic 



This would be most effective in combination with observa- 
tions constraining very early optical flashes, which have been 



field, which is also amplified by the shock ( Rees & Meszaros 



1992 Meszaros & Rees 19971. Radio observations are cru- 



cial for pinning down the evolution of the peak flux and 
two of the three characteristic frequencies of the broad- 
band synchrotron spectrum: the peak frequency and the syn- 
chrotron self- absorption frequency (with the third one being 
the cooling- break frequency). The radio emission from the 
blast wave typically peaks much later than at optical and X- 
ray wavelengths, on a time-scale of days to weeks depending 
on the observing frequency (Chandra & Prail||2012b , and in 



some cases can be observed months or even years after the 
initial GRB explosion (e.g. Frail et al.|2000a van der Horst 



et al.|2008 >. Although radio afterglows are quite faint, with 
the majority having sub-millijansky peak fluxes ( Chandra! 



& Frail 20121, a large number of these have well-sampled 
lightcurves at timescales from days to years (see Figure [I]). 

There are a few GRBs for which the radio observations 
started within the first day after the gamma-ray trigger, and 
some resulted in the detection of an early peak (e.g. |Kulkarni| 
et al.|1999 Frail et al.|2000b I . The favoured explanation for 
these early radio flares is emission from the 'reverse shock,' 
a phenomenon arising from dynamics at the shock front. 
As the shell of shocked matter slows and the surrounding 
medium begins to influence the dynamics, the interaction 
region becomes bounded by two shock fronts: one propagat- 
ing into the unshocked external medium ('forward shock'), 
and a second, 'reverse shock' propagating backwards into 
the relativistic flow (as seen from the flow rest frame; see 
e.g. |Kobayashi et al.|1999 l. This model may be complicated 
further if there are pre-existing shock fronts in the external 
medium due to stellar winds (Pe'er & Wijcrs 2006). Emis- 
sion from the reverse shock region may give rise to high flux 
levels across the frequency spectrum at much earlier times 
than the typical afterglow, and it has been predicted that 
radio flares due to the reverse shock are more likely to occur 
within the first hours after the GRB ( jMelandri et al.|2010| ). 

Although a reverse shock should in principle be formed, 
the brightness may vary considerably, perhaps resulting in 
only a small fraction of GRBs displaying radio flares. How- 
ever, this phenomenon is poorly constrained due to the 
scarcity of radio observations soon after the prompt gamma- 
ray emission. With the programme at AMI discussed in this 
paper we will be able to systematically probe this part of ob- 
servational parameter space and look for radio flares caused 
by reverse shock emission. The detections of, or upper lim- 
its on, early peaks at radio frequencies are important for 
putting constraints on the physical parameters of the blast 
wave, for instance its energy, the energy in electrons, and 
the magnetic field. Furthermore, they can help solve an 
outstanding issue in GRB jet formation, namely whether 
the jet is Poynting flux dominated or a baryonic outflow, 
by constraining the energy emitted as the relativistic flow 
begins to interact with the external medium (Piran 2000). 



searched for systematically for some years now (see Routing 
et al.|2006 and references therein) 



Besides the forward shock and reverse shock in GRB 
jets, there have been predictions of alternative sources of ra- 
dio emission in the case of short GRBs (those with prompt 
gamma-ray emission durations of less than ~2 seconds, |Kou-| 
veliotou et al. 19931. One favoured model for short GRBs 



is the merger of a binary system consisting of two neu- 
tron stars or a neutron star and a black hole. These merg- 
ers are thought to produce mildly to sub-relativistic out- 
flows, resulting in radio emission on timescales of weeks to 
years ( |Nakar fc Piran||2011| |Piran et al.||2012[ ). On shorter 
timescales the predictions are more uncertain, one possi- 
ble outcome being the generation of pulsar-like coherent 
emission by the rapidly rotating magnetar that is produced 
shortly after the binary merger ( |Pshirkov fc Postnov|20"To| ). 

There may even be unpredicted mechanisms that pro- 
duce early radio emission, and this follow-up programme 
will provide a systematic search for these at ~ 15 GHz fre- 
quencies. Intriguingly, recent observations at 1.4 GHz made 
with a 12 m single dish antenna at the Parkes radio obser- 
vatory have suggested the presence of short coherent bursts 
of radio emission within the first 20 minutes after the GRB 
onset, hinting at the presence of just such radio phenomena 
occurring at early times ( |Bannister et al.|2012| |. 

2.2 Sample sizes and selection bias 

The sample of GRB radio frequency observations presented 
in Chandra & Frail ( 2012 1 — albeit large — is far from com- 



plete, being both strongly sensitivity limited and biased due 
to target selection effects. The observations that have been 
made, and are still being carried out, at different facilities 
(e.g. Jansky Very Large Array (JVLA), Westerbork Syn- 
thesis Radio Telescope, Australia Telescope Compact Array, 
Giant Metrewave Radio Telescope) are restricted to a subset 
of possible GRB triggers due to limited available observing 
time. Subsequently only a small fraction of GRBs have been 
followed up, mostly those that have already shown interest- 
ing behaviour in the optical or X-ray bands. The same con- 
straints on observing time also restrict systematic follow-up 
of GRB non-detections, with a radio non-detection in the 
first few days after the burst often resulting in no further 
radio follow-up, despite the fact that the peak flux of the 
radio emission will often occur at later times. 

The recent upgrade of the JVLA allows for much deeper 
observations, which will help significantly in probing GRB 
radio afterglows at fainter flux levels. However, the num- 
ber of GRBs that the JVLA and other facilities will observe 
is still limited, as illustrated by the fact that of the eleven 
GRBs observed and reported upon here, only five have been 
followed up by other radio telescopes (and the observations 
made known to the follow-up community). Although AMI is 
not as sensitive as the JVLA, this programme will provide 
a uniformly selected sample which can be used for compar- 
ison with GRB follow-up at other wavelengths. With sys- 
tematic follow-up we should also detect many 15-GHz radio- 
afterglow peaks at later times when other observatories have 
stopped observing. 

In summary, after the first year(s) of this project, we 
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will have a uniquely large sample of early time radio obser- 
vations of GRBs, providing important new constraints on 
short (minutes to hours) timescale radio emission. We will 
also provide a systematically obtained and rapidly reduced 
dataset, flagging up interesting GRBs for deeper observation 
both with AMI and other observatories as appropriate. 



THE ARCMINUTE MICROKELVIN 
IMAGER LARGE ARRAY 



The Arcminute Microkelvin Imager Large Array (Zwart 
et al. 2008 I is a synthesis telescope composed of eight 



equatorially-mounted 12.7m dishes sited at the Mullard Ra- 
dio Astronomy Observatory at Lord's Bridge, Cambridge. 
The telescope observes in the band 12.0-17.9 GHz with eight 
0.75 GHz bandwidth channels. In practice, the two lowest 
frequency channels are not generally used due to a lower 
response in this frequency range and interference from geo- 
stationary satellites, which leaves an effective bandwidth of 
13.5-17.9 GHz and a central frequency of 15.75 GHz using 
channels 3-8. After calibration, the phase is generally stable 
to 5° for channels 4-7 and 10° for channels 3 and 8. The 
telescope measures a single linear polarisation (I + Q) and 
has a flux density sensitivity of ~ 3mJy s _1//2 . Flux density 
calibration follows the Perley-Butler 2010 scale (R. Perley, 
private communication). AMI-LA data reduction has been 



described extensively in previous works (see e.g. Zwart et al. 
|2008[ ) and the flux calibration is accurate to better than 5 per 
cent ( |Scaife et al.|2008| |Hurley- Walker et~al]|2009[ ). During 
normal operation, the AMI-LA performs continuous execu- 
tion of observations scheduled as entries in a queue, with 
observing programmes often set up for several days in ad- 
vance. Implementation of rapid response observing for this 
GRB study, as part of a new target-of-opportunity capability 
for AMI, involves automatic insertion of observing requests 
into the queue, if necessary displacing existing entries or an 
observation in progress. 



4 AUTOMATED RESPONSE SYSTEM AND 
OBSERVING POLICY 

The problem of automated transient detection and follow- 
up is multifaceted, with potential for detailed focus in many 
subfields such as detection, classification, resource alloca- 
tion, etc. For the purposes of this paper we restrict our 
discussion to three broadly defined sub-systems: the dis- 
tribution network, observation request triggering, and the 
automated telescope response. 



4.1 Transient alerts distribution via the VO Event 
network 

NASA uses the Gamma-ray Coordinates Network (GCN, 
Barthelmy et al. 1998 1 to rapidly distribute information 



about GRBs from satellites such as INTEGRAL, Swift and 
Fermi. Until recently, the most common medium was a short 
piece of plain text conforming to a specified format, dis- 
tributed via a dedicated central server / client network. 
While fast and effective for this purpose, this approach will 
not encompass the range of hierarchical data structures that 



are envisioned for describing the wealth of transient events 
detected by a new generation of astronomical facilities such 
as LOFAR, Pan-STARRS, GAIA, LSST etc. A robust net- 
work of many facilities will also require a more flexible and 
adaptable distribution model, i.e. one that is largely decen- 
tralized. 

The VOEvent standard (see footnot^] |Williams et al 



2006 1 has been developed with these needs in mind, pro- 



viding a flexible data format specified as an XML schema. 
The standard is transport independent (allowing for multi- 
ple distribution methods), allows for a distributed network 
of clients and servers, and has been previously employed 
by projects such as eStar, the Catalina Real Time Survey, 
and Sky Al ert ( |Allan et al.|2006||Drake et al.|2009||Williams 
et al.|2009[ ). GCN has now started distributing VOEvents (S. 
D. Barthelmy , private communication) using the VOEvent 
transport protocoj^Jwhich is an interoperable standard for 
VOEvent distribution. 

We tested two freely available software packages im- 
plementing the VOEvent transfer protocol: 'Dakota^ind 
'Cornet]^] Both provide means of generating, distributing 
and subscribing to streams of VOEvents. For this work, 
we use Comet, primarily because it is implemented in the 
Python programming language. The packages are dedicated 
to efficient and robust networking; processing of the infor- 
mation enclosed in the packets is delegated to an external 
program of the user's choice. In the current system imple- 
mentation we monitor for new GRB events via a direct con- 
nection to the GCN VOEvents service, and pipe received 
VOEvents to a custom software package described below. 



4.2 VOEvent monitoring and trigger generation 

The VOEvent identification and response routines were im- 
plemented using a Python package written to facilitate such 
tasks: 'pysovo,' which is available under an open-source li- 
cence from https://github.com/timstaley/pysovo. The 
package makes it easy for the user to extract information 
from the VOEvent XML packets; implement decision logic to 
determine whether an observation request will be triggered; 
determine whether the source is currently observable from 
a given observatory; generate requests formatted to custom 
observatory specifications, and send notification emails and 
SMS (i.e. text to mobile phone) alerts to designated recip- 
ients. The package is designed to be easily extensible, in 
order to allow interaction with multiple sources of transient 
astronomical alerts and multiple follow-up facilities. 

For this programme, we perform VOEvent monitor- 
ing and trigger generation using instances of Comet and 
pysovo running on a dedicated machine in Southampton. 
Each VOEvent has an ID string, referred to as an 'Inter- 
national Virtual Observatory Resource Name,' (IVORN), 
which uniquely identifies a VOEvent in much the same way 
that a URL identifies a resource on the internet. We use a 



1 'Sky Event Reporting Metadata Version 2.0,' Rob Seaman et 
al, 2011: http://www.ivoa.net/Documents/VOEvent/index.html 

2 'VOEvent Transport Protocol,' A. Allan and K.B. Denny: 
http: //www. ivoa.net/Documents/Notes/VOEventTransport/ 

° 'Dakota,' R.B. Denny: http://voevent.dc3.com 
4 'Comet,' J.D. Swinbank: http://comet.transientskp.org 
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simple filter based on the IVORN to identify the VOEvent 
information packets pertaining to new Swift GRB alerts. 
Once relevant VOEvents have been identified, the coordi- 
nates are extracted, checked against our observation policy, 
and used to generate a specially formatted observation re- 
quest email which is then sent to an address at the Mullard 
Radio Astronomy Observatory, Cambridge. 

4.3 Automation of the AMI telescope 

At the telescope level, incoming request emails to a ded- 
icated account are picked out by an email filter based on 
the subject line and a template contained in the message 
body. The request email is then passed to a custom program 
written to respond to AMI target of opportunity requests, 
'rqcheck,' running on the AMI- LA control computer. The 
email is also forwarded to observatory members, who can 
act as request moderators if necessary. 

The 'rqcheck' process validates the email, checking syn- 
tax and keyword values, and then proceeds to: 

• Check whether target is currently on-sky and not ob- 
scured by the Sun or Moon. Sets start and stop times for 
the next possible timing within the next 24 hrs, according 
to the timing specified in the request (ASAP, next transit, 
or specified sidereal time). 

• Find the nearest phase calibrator from VLBA^] and 
JVAS ( |Patnaik et al.|1992[ ) catalogues. 

• Check the availability of the telescope, e.g. in case of 
priority observing or engineering work. 

• Construct an entry for the observing queue. This will 
always be a standard pointing observation with interleaved 
calibrator, using the default integration times for either ar- 
ray. 

• Insert this entry in the queue, displacing any previously 
scheduled observations if necessary. 

• Append the request to a log file. 

• Notify the requester of the outcome by reply e-mail. 



4.4 Trigger and follow-up policy 

For observations with AMI we employ a naive trigger policy, 
whereby all Swift GRB alerts with a declination above -10° 
cause an observation request. This simple strategy is optimal 
in this case as we have access to sufficient observing time 
to allow systematic follow-up of all Swift GRB triggers in 
the northern hemisphere. We are employing a logarithmic 
follow-up schedule for all triggers, observing promptly after 
the initial burst, after 1-2 days, 3-4 days, 1 week, 2 weeks, 
and finally at 1 month. Follow up of detected sources will 
continue to be evaluated manually. Our standard integration 
time for each observation is currently one hour. 



5 DATA REDUCTION AND ANALYSIS 

When undertaking systematic observations such as those en- 
tailed by a GRB follow-up programme, it is desirable to have 
a fully automated reduction and analysis pipeline. With this 



goal in mind we have developed a procedure requiring min- 
imal human invervention, which we plan to fully automate 
in the near future. 



5.1 Data transfer and image synthesis 

After observations have been made, the raw data are trans- 
ferred from Cambridge to Southampton for further analy- 
sis. The data are then reduced using a procedure scripted in 
Python which performs the following tasks: 

• Identifies set of observations pertaining to each GRB. 

• Runs the AMI software tool 'reduce' on each obser- 
vation to perform flagging for interference, shadowing and 
hardware errors, apply gain calibrations and synthesise the 
frequency channels to produce calibrated visibility data in 
m^FITS format. 

• Extracts certain metadata which can be used to flag 
observations of poor quality (e.g. system temperature esti- 
mate, RFI flagging percentages). 

• Generates images using the casapy 'clean' algorithm, 
including a 'dirty' map and PSF map, for visual inspection. 



5.2 Image time-series analysis 

Once a time-series of images has been synthesised these are 
visually checked for for data quality, artefacts, and sources 
coincident with the known GRB position. We then attempt 
to extract a lightcurve by running custom source- fitting soft- 
ware developed for LOFAR data analysis (Sprccuw 2010[ ). 
This allows for source fitting at a known location, in addition 
to 'blind' source finding and extraction above a given thresh- 
old. To ensure that marginal detections are not missed, we 
use the source-fitting software to attempt an extraction at 
the location given by the enhanced Swift X-ray Telescope 
reduction (Burrows et al. 20051; except where this is un- 



available (e.g. in the case of GRB120403A) and we resort to 
the Swift-BAT refined position. In the case of a detection, 
the code attempts to fit an elliptical Gaussian profile, and 
errors on the fit are calculated according to the formulae of 



5 NRAO VLBA catalogue: |http : //www ■ vlba . nrao~ edu/ astro/ 
calib/ 



Condon (19971. When no detection is made, we assign an 
upper limit of three times the image RMS level. 



6 RESULTS 

6.1 Response times 

We can estimate system response times by inspecting times- 
tamps on the GCN notices and the AMI observation request 
logs. The typical response time for a target which is im- 
mediately observable with AMI is approximately 5 minutes 
from Swift-BAT trigger time to AMI taking data, and we 
can break this down into contributing factors. First, Swift 
GRB positions are broadcast from GCN, with a typical time 
delay of 7-30 seconds after the trigger time (approximately 
0.1 seconds of which is due to GCN — S. Barthelmy, private 
communication). These are received by our VOEvent node 
in Southampton, processed, and an email is sent to the AMI 
control system. Round trip time from NASA-GCN to AMI 
is consistently 4-6 seconds, including the time required to 
authenticate with the email server. Alerts distributed solely 
via VOEvent transfer protocol will be relayed with an even 
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GRB ID 


Hours 
since burst 


3<t upper limit 
(mJy) 


\j JrL £3 ± A U o U O J\. 


U.U / 


(1 QIC 


GRB120308A 


0.08 


0.164 


GRB120311A 


0.07 


0.230 


GRB120320A 


14.32 


0.196 


GRB120324A 


0.07 


0.218 


GRB120326A 


7.37 


0.430 


GRB120403A 


7.48 


0.238 


GRB120404A 


11.32 


0.223 


GRB120422A 


6.41 


0.616 


GRB120514A 


50.91 


0.211 


GRB120521C 


0.24 


0.302 



Table 1. Delay between GRB trigger timestamp and the start 
time of our initial observation, for each GRB followed up. No ra- 
dio sources at the position of the GRB were detected at the 3<r 
level in our preliminary analysis of these data, and we list corre- 
sponding upper limits for the flux density averaged over the one 
hour observation period. Note that GRB 120326 A has a marginal 
detection at 0.337 mJy in the initial observation, ~ 7 hours after 
the trigger timestamp. 



shorter delay. If the target is on the sky, AMI then begins 
to slew almost immediately. The maximum slewing speed of 
the telescope is approximately 15° per minute, which means 
that response times for immediate observation will be of the 
order of a few minutes; currently a 4-minute period is allot- 
ted to ensure that the telescope is on source before taking 
data. 

Immediate observations with AMI are not always pos- 
sible — the telescope may be undergoing maintenance, ex- 
periencing unfavourable weather conditions, or the source 
may simply be over the horizon. However, even initial ob- 
servations delayed by a few hours are faster than the large 
majority of previously catalogued radio follow-up, as illus- 
trated in Figure [l] 



6.2 Radio follow-up lightcurves 

Delays between each GRB trigger and our initial follow-up 
observation, and corresponding 3a upper limits on the flux 
density, are given in Table [l] We detect a radio afterglow 
with confidence level above 5a for GRB 120326A, as shown 
in Figure [2] The remaining ten have lightcurves consistent 
with non-detection. Specifically, in this preliminary analysis 
we do not detect any radio emission from the five GRBs 
observed within the first six hours post-GRB trigger. This 
places upper limits of less than 1 mJy on the 15-GHz flux 
density averaged over the one hour observation period, but 
only weakly constrains possible emissions on much shorter 
timescales. More detailed analysis of the data and upper 
limits is planned for future work, as discussed in Section [7 
A full listing of all observations is given in Appendix \K 
The BAT and BAT-XRT lightcurves shown in Figure[2]were 
obtained from the automated burst analyser page at the UK 
Swift Science Data Centre Website ( Evans et al.|2010| using 
4er significance bins and plotting the observed flux density 
at 10 keV. 
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Figure 1. GRB radio follow-up delay histograms. Top: A his- 
togram depicting radio observations of GRBs made with the VLA. 
Light grey bars show counts for all observations, while medium 
grey bars (overlaid) show counts for the first observation of any 
given source. Bottom: A close up depicting the observation counts 
for initial radio observations of a GRB source undertaken less than 
16 hours after the burst. Medium grey bars depict VLA observa- 
tions. Dark grey bars, overlaid, depict observations presented in 
this paper (there are no catalogued VLA observations within the 
30 minute period immediately after a burst). Data on VLA ob- 
servations is reproduced from |Chandra fc Frail] ( |20 12[ ) . 



7 DISCUSSION AND FUTURE WORK 

In summary, this paper presents the first results from our 
programme of GRB follow-up observations using AMI-LA, 
demonstrating the success of the fully automated rapid- 
response follow-up system we have developed. Future obser- 
vations will provide the first large sample of early-time radio 
observations of GRBs, placing tight constraints on theoret- 
ical models, while our monitoring strategy will result in a 
sample of radio afterglows unbiased by selection effects. 

There are a number of enhancements to the AMI-LA 
response system we hope to implement in the near future. 
These include: 

• Notifying the wider astronomical community when 
AMI-LA follow-up observations have been made, both 
through the NASA-GCN and the VOEvent network. 
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BAT-XRT data of 1 20326A 
BAT: Black — XRT: WT: Blue; PC: Red 
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Figure 2. Lightcurves for GRB120326A. Left: Swift BAT and XRT measurements of the gamma- and X-ray emission. Right: Radio- 
afterglow lightcurve. Star points with errorbars represent 15-GHz observations made with AMI as part of this programme. Circle points 
represent 3<r blind detection thresholds — a connecting dashed line is overplotted to guide the eye. The diamond point at 5.45 days 
represents a J VLA detection at 22 GHz reported via GCN circulars | |Laskar et al.||2012| l . Vertical dashed line in both plots represents 
the time at which AMI first started taking data on the GRB source position. Note that GRB 120326A was observed after a much longer 
delay than is possible with the automated response system (~7 hours rather than the typical 5 minutes). 



• Continued refinement of the data reduction process, in- 
cluding automated data transfer and initiation of the data 
reduction process. We will also implement time-slicing and 
analysis of the data on timescales shorter than one hour, 
to look for short-duration, high intensity emission. This is 
a particularly exciting avenue of investigation given recent 
results reported in Bannister et al.| ( |2012[ ). 

• Automated discovery and user notification of GRB af- 
terglow candidates, to enable rapid human intervention and 
detailed follow-up. 

We note that a human-verified dataset of modest size such as 
this provides an excellent test bed for automated transient- 
detection and classification tools. We are working closely 
with the LOFAR Transients Key Project team on software 
to this end ( |Fender et aT1|2006| [Swinbank et"al~|2011| >, and 
intend to employ it to automate much of the data curation 
currently done manually in the course of these follow-up 
observations. 

This work was made possible in part through ERC fund- 
ing for the 4 Pi Sky project]^] In a wider context, this work 
provides a first step toward the longer term goals of that 
project, such as automated data reduction and distribution 
of machine readable data products (not just notifications) to 
the astronomical transients community, and use of such data 
products to enable rapid, automated decisions on allocation 
of telescope resources for transient follow-up. 

It is hoped that the software described in this paper 
will be useful to the wider astronomical transients commu- 
nity in working toward these aims. We also plan to imple- 
ment follow-up protocols in collaboration with other obser- 
vatories. Notably, we are collaborating with the LOFAR-UK 
station ( Fender|2010 | which will provide extremely rapid re- 
sponse due to the way in which LOFAR stations are made up 

6 http://www.4pisky.soton.ac.uk 



of phased arrays (and hence have no moving parts). We plan 
to start taking data with LOFAR-UK in the near future. 
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APPENDIX A: TABLE OF OBSERVATIONS 

A full list of all AMI-LA 15-GHz observations presented in 
this paper is given in table [AT] 
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Flux 


Image 


Start 


Days 


den- 


std. 


date 


since 


sity 


dev. 


(UTC) 


burst 


(mJy) 


(mJy) 


GRB120305A 




12/03/05 19:41:57 


<0.01 


<0.316 


0.105 


12/03/06 14:20:54 


0.78 


<0.158 


0.053 


12/03/07 15:26:46 


1.83 


<0.268 


0.089 


12/03/08 15:52:45 


2.84 


<0.171 


0.057 


12/03/09 18:08:27 


3.94 


<0.227 


0.076 


12/03/10 16:14:49 


4.86 


<0.193 


0.064 


12/03/11 17:50:36 


5.93 


<0.293 


0.098 


12/03/12 14:57:08 


6.81 


<0.272 


u.uyi 


12/03/13 17:17:49 


7.90 


<0.218 


0.073 


12/03/15 17:44:51 


9.92 


<0.226 


0.075 


12/03/16 18:00:52 


10.93 


<0.163 


0.054 


12/03/19 17:29:07 


13.91 


<0.203 


0.068 


12/03/22 13:08:01 


16.73 


<0.221 


0.074 


12/03/27 16:17:47 


21.86 


<0.197 


0.066 


12/04/01 16:18:04 


26.86 


<0.194 


0.065 


12/04/05 16:42:14 


30.88 


<0.169 


0.056 


GRB120308A 




12/03/08 06:18:20 


<0.01 


<0.164 


\J.VJOO 


12/03/09 02:16:03 


0.83 


<0.144 


0.048 


12/03/13 03:20:06 


4.88 


<0.267 


0.089 


12/03/16 06:02:50 


7.99 


<0.183 


0.061 


12/03/19 05:11:09 


10.96 


<0.199 


0.066 


12/03/22 03:54:32 


13.90 


<0.217 


0.072 


12/03/25 03:27:46 


16.88 


<0.210 


0.070 


12/04/02 03:46:11 


24.90 


<0.169 


0.056 


12/04/11 02:25:55 


33.84 


<0.190 


0.063 


GRB120311A 




12/03/11 05:37:37 


<0.01 


<0.230 


0.077 


12/04/05 04:39:13 


24.96 


<0.186 


0.062 


12/04/08 03:37:33 


27.92 


<0.180 


0.060 


GRB120320A 




12/03/21 02:15:44 


0.60 


<0.196 


0.065 


12/04/05 01:24:45 


15.56 


<0.212 


0.071 


GRB120324A 




12/03/24 06:03:17 


<0.01 


<0.218 


0.073 


12/04/05 05:49:01 


11.99 


<0.175 


0.058 


12/04/08 04:47:22 


14.95 


<0.240 


0.080 


12/04/16 08:25:13 


23.10 


<0.184 


0.061 


12/04/23 03:58:21 


29.92 


<0.329 


0.110 







Flux 


Image 


Start 


Days 


den- 


std. 


date 


since 


sity 


dev. 


(UTC) 


burst 


(mJy) 


(mJy) 


GRB120326A 




12/03/26 08:42:58 


0.31 


0.337 


U.14o 


12/04/02 04:55:59 


7.15 


0.771 


U.Uo4 


12/04/04 07:37:40 


9.26 


0.666 


U. loy 


12/04/08 07:31:55 


13.26 


0.757 


\J.Zl ( 


12/04/16 07:15:25 


21.25 


0.626 


0.093 


12/04/28 07:03:08 


33.24 


0.495 


0.161 


GRB120403A 




12/04/03 08:34:27 


0.31 


<0.238 


U.U ( y 


12/04/04 16:06:16 


1.63 


<0.205 


U.UOo 


12/04/06 16:08:23 


3.63 


<0.204 


0.068 


12/04/10 14:52:49 


7.57 


<0.212 


0.071 


12/04/16 15:34:03 


13.60 


<0.254 


0.085 


12/05/01 14:50:02 


28.57 


<0.349 


0.116 


GRB120404A 




12/04/05 00:09:57 


0.47 


<0.223 


U.U / 4 


12/04/08 00:18:06 


3.48 


<0.237 


n ri70 
u.u / y 


12/04/11 00:11:17 


6.47 


<0.169 


0.056 


12/04/17 01:57:21 


12.55 


<0.214 


0.071 


GRB120422A 




12/04/22 13:36:43 


0.27 


<0.616 


U.zUO 


12/04/24 16:17:23 


2.38 


<0.465 


U. loo 


12/04/26 16:19:30 


4.38 


<0.368 


U. IZo 


12/04/29 17:17:31 


7.42 


<0.239 


0.080 


12/05/05 17:53:45 


13.45 


<0.228 


0.076 


12/05/20 21:19:03 


28.59 


<0.456 


0.152 


GRB120514A 




12/05/16 04:07:39 


2.12 


<0.211 


U.U / u 


12/05/19 03:01:00 


5.08 


<0.460 


U. lOo 


12/05/22 01:54:22 


8.03 


<0.279 


u.uyo 


12/06/03 23:58:25 


20.95 


<0.619 


0.206 


12/06/11 03:20:16 


28.09 


<1.488 


0.496 


GRB120521C 




12/05/21 23:36:44 


0.01 


<0.302 


0.101 


12/05/23 21:02:16 


1.90 


<0.793 


0.264 


12/05/25 20:54:25 


3.90 


<0.307 


0.102 


12/05/28 21:02:34 


6.90 


<0.328 


0.109 


12/06/04 18:30:23 


13.80 


<0.216 


0.072 


12/06/23 18:50:25 


32.81 


<0.342 


0.114 



Table Al. 15-GHz measurements of the first eleven GRB sources observed with AMI-LA. For sources 
detected (i.e. all except GRB 120326A), the flux density is given as an upper limit at the 3a level. 



rhere radio emission was not 
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